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Abstract 

Taxonomic investigation of spined loaches from Dalmatia and Herzegovina was conducted on specimens from 14 localities. 
The results of the detailed morphological investigations were combined with genetic data (based on one mitochondrial and 
two nuclear genes) in order to resolve the taxonomic status of each Cobitis population. Among the investigated features of 
external morphology, the appearance of spots on the caudal fin base turned out to have the greatest diagnostic value. 
Furthermore, the number of branched fin rays enabled the discrimination of several species. No morphometric character 
alone could ensure determination of any Cobitis species. Nevertheless, groups of populations that are more similar in their 
body shapes correspond to mitochondrial phylogenetic lineages. Based on molecular genetic markers, Dalmatian and 
Herzegovinian spined loaches form independent lineages inside the Adriatic phylogenetic group. Mitochondrial DNA 
phylogenetic reconstruction revealed six monophyletic lineages, corresponding to six species distributed in the investigated 
area. The population distributed in IVlostarsko blato karstic field in Bosnia and Herzegovina is described as a new species 
based on a unique combination of morphological characters: a single triangular Canestrini scale; usually 5V2 branched anal 
fin rays, 6V2 branched dorsal fin rays, 14 branched caudal fin rays; no spots in the surface pigmentation layer on the caudal 
fin base; scales on the body very small. 
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Introduction 

Recognition and delimitation of species, the main objectives of 
systematics, are of crucial importance for implementation of 
adequate strategies for biodiversity conservation. Exclusive 
reliance on morphological characters, that were traditionally used 
to identify species, turned out to underestimate diversity by failing 
to detect cryptic taxa [1]. Genetic data are frequently used to 
delimit species, when conclusions about species status are based on 
an exclusivity criterion, such as reciprocal monophyly or degree of 
genetic clustering [2], [3]. However, those methods are also 
subjective and fail to account for several phenomena that cause 
gene tree incongruence with species tree [1]. Finally, recent 
implementations of coalescent model-based approaches allow 
conduction of probabilistic tests of species limits [4], [5]. Although 
some authors consider that integrative methods will not yield a 
clear result in cases when cryptic species are present [1], the 
critical review of the results of all mentioned methods may have 
the greatest resolution power in delimiting species, especially 
in cases when phenomena such as hybridization are more 
widespread, species are morphologically very similar and/or 



divergence happened more recently. Among such, "more difficult" 
taxonomic cases are freshwater spined loaches. 

Even though allopatric speciation has been highlighted as 
responsible for divergence of many freshwater species, there are 
evidences of sympatric/parapatric speciation modes [6], [7], as 
well as contribution of genetic mechanisms other than selection 
(i.e. drift and founder events) [6] in freshwater fish studies. Fishes 
in general present a very diverse array of speciation scenarios 
(reviewed in [8]). 

The genus Cobitis (spined loaches) comprises about 60 species of 
freshwater fishes, [9], [10], [11] that are distributed throughout 
the temperate zone of Europe and Asia. Due to the morphological 
similarity of different Cobitis species, and the existence of sibling 
species [12], their taxonomic dilferentiation is complicated and 
their systematics not completely resolved. Sibling species are 
defined as a cryptic sister species; two species that are the closest 
relatives and have not been distinguished from one another 
taxonomicaUy [13]. Recent descriptions of new species and the 
discovery of significant diversity among spined loaches has been 
made possible by a combination of detailed morphological 
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investigations (e.g. [1 1], [12], [14], [15]), karyological analyses (e.g. 
[16], [17]) and analysis of DNA markers (e.g. [18-22]). 
Dalmatia, which includes rivers in the Adriatic watershed in 

Croatia and Bosnia and Herzegovina, is well known for its high 
degree of endemism of freshwater fishes [23-25]. Several 
molecular genetic investigations of spined loaches in this region 
[18], [21], [22] have revealed that it is inhabited by Cobitis 
populations belonging to the Adriatic phylogenetic group, a 
distinct monophyletic group. It is currently thought that the 
Adriatic watershed in Croatia is inhabited by five different Cobitis 
species (Figure 1): C. jadovaensis Mustafic & Mrakovcic, 2008, a 
recently described species living only in the small karstic river 
Jadova [11]; C. bilineata Canestrini, 1865, which inhabits the 
Zrmanja River in Croatia [21], [22], but is also distributed in 
Slovenia, Italy, France, Switzerland and Spain [9]; C. dalmatina 
Karaman, 1928, endemic to the Cetina River in Croatia [26], 
[14]; C. illyrica Freyhof & Stelbrink, 2007, another recentiy 
described species with its only known locality being Prolosko blato, 
a small lake in the Imotsko Polje karstic field [15]; and C. narentana 
Karaman, 1928, supposedly distributed in basin of the Neretva 
River, namely in the lower courses and tributaries of the Neretva 
River, as well as in the Bacinska Lakes and Matica River [27]. 
Recent descriptions of two species from southern Croatia (C. 
jadovaensis and C. illyrica) were based on morphological characters 
[1 1], [15]. In the Adriatic watershed in Bosnia and Herzegovina 
spined loaches have been recorded in the Neretva basin (including 
Hutovo blato wetland and Mostarsko blato karstic field), in 
Trebisnjica River, and also in Krenica Lake, which is located in 
Bekijsko karstic field [28] (Figure 1). Whereas Sanda et al. [28] 
suggested the existence of another Cobitis species in addition to C. 
narentana in this area of Bosnia and Herzegovina, neither the 
taxonomic status of the populations nor the distribution of each 
species have been ascertained. 

Although previous investigations recognized the region of 
Dalmatia and Herzegovina as an area inhabited by a large 
number of Cobitis species comprising many endemics, the diversity 
of spined loaches in the area is probably still underestimated, while 
data on intraspecific diversity being especially scarce. The greatest 
problem with the investigations performed so far is that all of them 
were based either on morphological characters or on DNA 
markers, and none on both, and usually on ver^' few specimens 
from one or a few populations. This is especially true for the 
description of C. illyrita [15] where the sampk; size analyzed is too 
low to represent the morphological diversity and distribution of the 
species. Original descriptions of C. dalmatina and C. narentana [26] 
are also problematic, consisting of a few very general sentences in 
the original text. Other investigations gathered under one 
supposed taxon specimens that, based on today's knowledge, 
belong to more than one species. For example, in the investiga- 
tions of Schneider et al. [14], [27] in the sample of C. narentana, C. 
illyrica specimens were also included. 

The aim of the present investigation was to determine the 
taxonomic status, phylogenetic relationships, level of intraspe- 
cific and intrapopulational diversity and the distribution of 
spined loaches in the region of Dalmatia and Herzegovina. 
Furthermore, based on large samples from more than one 
population (where possible) for each species, we tested the 
diagnostic markers traditionally used in Cobitis species determi- 
nation and ascertained the diagnostic value of each morpho- 
logical character. 



Materials and Methods 

Ethics Statement 

This investigation was conducted entirely in accordance with 
ethical standards and Croatian legislation. The work was approved 
by the Ethical Committee of the Faculty of Science, University of 
Zagreb. The permission for sampling on all localities in Croatia 
was issued by the Ministry of Agriculture, Forestry and Water 
Management of the Croatia. The permission for sampling on 
Hutovo blato, Trebisnjica, Krenica and Mostarsko blato was 
issued by the Ministry of Agriculture of Bosnia and Herzegovina. 
The sampling was conducted by electrofishing. The fishes were 
over-anaesthetized with tricaine-methanesulphonate (MS 222) and 
preserved in ethanol. All specimens are deposited in institutional 
collections (see List of examined material). 

Specimens were collected from 14 localities in Dalmatia and 
Herzegovina (List of examined material. Figure 1). All rivers and 
lakes where spined loaches have been recorded inside the 
investigated area were included in the study. Where possible, 
sampling was conducted at more than one locality in each river. 

In order to ascertain the taxonomic position and phylogenetic 
relationships of the Adriatic spined loaches, two t^-jx^s of analyses 
were employed - the investigation of morphology and the analyses 
of molecular phylogeny. Morphological analyses comprised the 
analyses of morphometry, meristics and external morphology. 

A total of 24 morphometric characters were measured using an 
electronic caliper to the nearest 0.01 mm: total length (TL), 
standard length(SL), lateral head length (c), preanus distance (Pan), 
preanal distance (aA), prepelvic (preventral) distance (aV), 
prepectoral distance (aP), predorsal distance (aD), caudal peduncle 
length (Ipc), length of dorsal (ID), anal (lA), caudal (IC), pectoral (IP) 
and pelvic (ventral) (IV) fins, distance between pelvic fins and anal 
aperture (Van), head depth (he), maximum body depth (H) and 
caudal peduncle depth (or minimum body depth, h), head width 
(laco), maximum body width (lac), distance between eyes (io), eye 
diameter (o), preorbital distance (prO), postorbital distance (poO). 
All measurements were made point to point by one author. Size- 
dependent variation was removed by allometric transformation of 
morphometric measures using the formula [29], [30]. 

Madj = M{SmSLof , 

where M^dj is the adjusted size-independent measurement, M the 
original morphometric measurement, SL^ the overall mean of the 
standard length for all fishes from each population, SL„ the 
standard length of the fish, and b the slope of the regression of log 
M on log SLj, based on all specimens from each population. 
Correlation coelHcients between transformed, adjusted measure- 
ments and standard length were calculated to confirm that the 
effect of size was removed. 

Morphometric ratios, i.e. percentage ratios of morphometric 
characters in relation to SL, c and H, as well as h/lpc were also 
calculated in order to enable comparison with data in the 
literature, especially in cases when certain morphometric ratios 
were used to separate species. 

To compare morphometric features between males and females. 
Student's t-test was employed. Statistical comparison among 
populations was conducted using analysis of variance (ANOVA) 
and principal component analysis (PCA), and was based on size- 
independent measurements. Software package STATISTICA 7.1 
was used for data analysis. 

The meristic characters assessed included the number of 
unbranched and branched fm rays in the dorsal, anal, pelvic, 
pectoral and caudal fins. The last two branched rays in dorsal and 
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Figure 1. Map of Investigated area with sampling localities mariced, as well as distribution of Cobitis species based on literature. 

Legend: /\ = C. jadovaensis; B = C. bilineata; C = C. dalmatina; D = C. illyrica; E = C narentana; 1 =Jadova; 2 = Zrmanja (HE Velebit accumulation); 
3 = Cetina (Blato); 4 = Neretva in Metkovic; 5 = Mislina; 6 = Norin; 7 = Modro ol<o; 8 = Bacinsl<a Lal<es; 9 = Matica; 10 = Prolosi<o blato; 11 =Trebisnjica; 
12 = Hutovo blato; 13 = Mostarsko blato; 14 = Krenica; A = Austria; H = Hungary; SLO = Slovenia; HR = Croatia; I = Italy; BIN = Bosnia and Herzegovina; 
RO = Romania; SRB = Serbia; MK = The former Yugoslav Republic of Macedonia. 
doi:1 0.1 371/journal.pone.0099833.g001 



anal fins, which articulate on a single pterygiophore, were scored 
as "1 V2". 

The external morphology was examined in detail in order to 
detect difTerences between populations and species and to reveal 
phylogenetic signals (defined as the degree to which phylogenetic 
relatedness among taxa is associated with their phenotypic 
similarity [31]) in their morphology. Attention was paid to 
characters traditionally used in Cobitis species difierentiation, 
including the appearance of spots on the caudal fin base (see 
[9], [11], [14], [15], [32]), the appearance of Gambetta zones on 
body sides (see [9], [12]), structure of the scales (see [12]) and the 
position and shape of the suborbital spine (see [14], [32]). 
Gambetta zones [33] are pigmentation zones comprised of various 
spots and blotches. Loaches of the genus Cobitis have four 
Gambetta zones located along the body [9], usually recorded as 
Z1-Z4, zone Zl being located dorsolaterally and Z4 ventrolater- 
ally. Saitoh & Aizawa [34] reported that zone Z4 consists of two 
pigmentation layers; surface layer located in the dermis and deeper 
layer in the horizontal septum and its surrounding tissue. Likewise, 
pigmentation located on the beginning of the caudal fin is also 
formed by two pigmentation layers [34]; surface layer in the 
dermal tissue and the deeper layer in the connective tissue 
surrounding the base of fin rays. In this investigation we 
differentiated two layers of pigmentation, unlike majority of 
authors reporting morphological characters of spined loaches (e.g. 
[II] [14] [15]). Thereafter, our results may not be completely 
comparable with others due to the lack of information whether 
pigmentation patterns reported by other authors refer to the 
surface layer only or to the combination of both layers. The 
suborbital spine (spinum suborbital^ is a modified ethmoidal bone, 
situated below the eye. 



An investigation of morphological features was conducted on 
226 fishes collected from 1 3 sites (List of examined material) and 
preserved in ethanol. However, only a few males were sampled so 
a normal distribution could not be confirmed at the following 
localities: Jadova, Neretva in Metkovic, Norin and Prolosko blato. 
Therefore, males from those locations were not included in the 
statistical comparative analyses. In Bacinska Lakes, despite 
intensive and widespread sampling and previous reports testifying 
to a large population of spined loaches (e.g. [27]), only two females 
were caught. Hence only molecular investigations were conducted 
on the specimens from these lakes. 

Molecular analysis was undertaken on 122 specimens collected 
from 14 localities (List of examined material. Table 1) and was 
based on three molecular markers: mitochondrial cytochrome b 
(cyt b), nuclear RAGl gene and the first intron of nuclear S7 gene. 

Total genomic DNA was extracted from fresh and deep-frozen 
fin tissue using a standard extraction product (DNeasy tissue kit, 
Qiagen). Polymerase chain reaction (PGR) amplifications were 
performed in a 50(xl reaction volume containing 25|J,1 of 
HotStarTaq Master Mix (Quiagen), 2|xl of each primer and 4^1 
of template DNA. Cyt b gene was amplified using primers LI 47 25 
and HI 6460 [35] and the following temperature regime: 15 min 
at 95°C; 35 cycles of 30 s at 94°C, 30 s at 50°C and 90 s at 72°C; 
7 min at 72"C. The primers and PGR protocols for amplification 
of RAGl and S7 genes are described in [22], [36] and [37]. 
Sequencing was carried out by Macrogen Service Centre (Seoul, 
South Korea) with internal primers H-COB_cyt638 and L- 
Cyp_425 [21] for cyt b. The same primers as for the PGR were 
used for sequencing of RAGl and S7 genes. 

Homologous regions of cyt b and RAGl genes were aligned 
manually against previously published sequences. Chromatograms 
and alignments were checked visually and were found to contain 
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no gaps or stop codons. Because of the presence of insertions and 
deletions in the S7 gene, those sequences were aligned with Clustal 
X [38]. Haplotype variants of nuclear genes in heterozygous 
individuals were reconstructed by a Bayesian statistical method 
implemented in PHASE 2.1 software [39], [40]. The phase 
reconstruction method (described in details in [39]) regards the 
unknown haplotypes as unobserved random quantities and aims to 
evaluate their conditional distribution in light of the genotype 
data, using Gibbs sampling, a type of Markov Chain Monte Carlo 
algorithm. The analyses were run five times for both nuclear genes 
with diffcrc'nt \alues of the seed of the pseudo-random number 
generator. Each run consisted of a burn-in-period of 100 followed 
by 1000 iterations. The gametic phase estimation was consistent 
through the runs according to goodness-of-fit measures. Nucleo- 
tide composition was examined by the homogeneity test of base 
frequencies implemented in PAUP (version 4.0bl0 [41]) for each 
data set, as well as the transition/ transversion ratio (ts/ tv). In order 
to test whether all mutations were selectively neutral, statistical 
tests D* and F* (proposed by Fu & Li [42]) and that of Tajima [43] 
were conducted using DnaSP v5 [44]. The same software was 
employed to estimate the recombination parameter, R [45] and 
the minimum number of recombination events, RM [46] , for each 
data set. For phylogenetic analyses, in addition to the sequences 
obtained in the present study, sequences available on GenBank 
were also included (Table 1). 

Pairwise comparisons of uncorrected sequence divergence (p- 
distances) of cyt b and RAGl were analyzed using MEGA version 
3.1 [47]. P-distances were not calculated for the S7 first intron 
because the presence of insertions and deletions of significant 
length (up to 1 1 base pairs) would lead to overestimates. 

In order to ascertain the phylogenetic position of Dalmatian and 
Herzegovinian spined loaches within the phylogenetic tree of 
European Cohitis, and to resolve the phylogenetic relationships 
among them, three different methods of phylogenetic reconstruc- 
tion were employed: maximum likelihood (ML), maximum 
parsimony (MP), both implemented in PAUP, and Bayesian 
inference (BAY), implemented in MrBayes (version 3.1.2 [48]). 
ML analysis was performed under the heuristic search option 
using the tree bisection-reconnection (TBR) branch swapping 
algorithm. For MP analysis, a heuristic s(;ar( h mode with 100 
replicates was used, with randomized input orders of taxa, and 
TBR branch swapping with all codon sites and nucleotide 
substitutions types weighted equally. Nonparametric bootstrapping 
(100 pseudo-replicates for ML, 1000 for MP, 10 additional 
sequence replicates) was used to assess branch support (BS). Each 
BAY analysis consisted of two simultaneous runs. For each, 
Markov Chain Monte Carlo was run four times for three million 
generations with trees sampled every 100 generations. The first 
20% of the sampled trees were discarded and Bayesian posterior 
probabilities (BPP) were estimated from the 50% majority-rule 
consensus tree of the retained trees. For ML and BAY analyses, 
the best-fitting model of molecular evolution was selected by 
hierarchical likelihood ratio tests using MODELTEST software 
(version 3.06 [49]). Phylogenetic analysis was performed on each 
investigated genetic marker independendy, due to differences in 
data sets composition (depending on the available sequences in 
GenBank) and differences in the mutational rate and phylogenetic 
performance of each investigated gene. Sequences of Sahanejewia 
romanica (Bagescu, 1943) [18], which belongs to the same family as 
Cobitis, and Cyprinus carpio L., 1 758, from a different family though 
the same order, were used as outgroups for cyt b phylogenetic 
reconstruction. For the RAGl data set, a s(^quence of 6". balcanka 
was used for rooting, while for S7, C. elongatoides and C. strumkae 
served as outgroups. 



In order to incorporate the possibility of horizontal gene transfer 
into the phylogenetic reconstruction of nuclear DNA, phylogenetic 
networks of nuclear haplotypes were constructed based on a 
median -joining (MJ) algorithm using Network 4.5.1.6. software 

(Fluxus Technology' Ltd.). 

To undertake a detailed phylogenetic reconstruction of spined 
loaches from Matica R., Prolosko blato, Bacinska lakes, Krenica 
and Mostarsko blato, and to confirm their taxonomic status, cyt b 
sequences from those localities were analyzed by a statistical 
parsimony method [50] under a 95'% connection limit, using the 
program TCS (version 1.21 [51]). 

As a complementation to already described morphological and 
phylogenetic analyses, we have also employed Bayes factors (BF) 
for comparison of phylogenetic hypotheses or models. That 
approach is one of recentiy described coalescent-based methods 
for statistical species delimitation using multi-locus DNA sequence 
data [5]. BF is calculated as the ratio of the marginal likelihoods of 
two models, which has the advantage of taking into account priors 
used in Bayesian analyses [52]. Even though reliance on Bayes 
factors as a species delimitation model tool has not been 
thoroughly explored [5], available examples approve their 
usefulness in model averaging [1], [5]. In this investigation we 
have used BF to test six phylogenetic models (scenarios) that were 
created based on the results of all other analyses. Model A 
corresponds to our conclusion that each mtDNA sublineage 
represents separate species (six species); in model B seven species 
were included because C. bilimata was considered as two species. 
Model C comprised six species, because C. illyrica and C. 
herzegoviniensis Buj & Sanda, sp. nov. were considered as one 
species; whereas for model D C. dalmatina and C. narentana are also 
considered as one species. Model E relies on the RAGl gene tree 
so it comprises four species (C.jadovaensis, C. narentana + C. bilineata, 
C. dalmatina, C. illyrica + C. herzegoviniensis Buj & Sanda, sp. nov.). 
Likewise, model F is based on the S7 first intron so it contains only 
three species [C. narentana + C. bilineata, C. dalmatina, C. Ulrica + C. 
herzegoviniensis Buj & Sanda, sp. nov. + C. jadovaensis). We used the 
coalescent-based species tree program *BEAST [53] and prepared 
multi-locus data set (comprising sequences of three genes) with 
three partitions. Analyses were performed under an uncorrelated 
lognormal relaxed molecular clock, whereas a Yule process was 
used for the species tree prior. Analyses were run for one million 
generations with 20% of the trees discarded as burn-in. Marginal 
likelihood scores, as well as Bayes factors were calculated in Tracer 
vl.5. [54]. Besides the harmonic mean estimator (HME), we have 
also employed path sampling (PS) and stepping-stone sampling 
(SS) for estimation of the marginal likelihood values of competing 
models. Namely, previous investigations demonstrated better 
performance of recently developed techniques PS and SS than 
HME by not overestimating the true marginal likelihoods [1] [5]. 
PS explores an almost continuous progression of distributions 
along a path from the posterior to the prior when calculating the 
marginal likelihood [55], whereas stepping-stone sampling bridges 
the posterior and prior distributions [52] [5]. For PS and SS 
calculations appropriate code [1] was added into XML files 
generated by BEAUti vl.7.4 [56]. 

Finally, nucleotide sites that represent fixed differences between 
Cobitis species (as proposed in this investigation), as well as shared 
polymorphisms, were detected using SITES [57]. 

List of Examined IVlaterial 

Institutional abbreviations: NMP, National Museum, Prague; 
PMF, Faculty of Science, Zagreb; PMR, Natural History 
Museum, Rijeka. 
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Table 1. Sampling localities, haplotype codes and GenBank accession numbers of specimens included in the phylogenetic 
analyses. 





locality 


species 


cyt b haplotypes 


RAG1 

haplotypes 


S7 haplotypes 


accession number 


reference 




Jadova 


C. jadovaensis 


JADl 


rJADl,2 


sJADl-4 


KJ487435, KJ487503, 
KJ487504, KJ487534- 
KJ487537 


this study 




Zrmanja 


C. bilineata 


ZRMl-5 


rZRMl-3 


sZRMl,2 


KJ487450, KJ487464, 
KJ487467, KJ487468, 
KJ487471, KJ487517- 
KJ487519, KJ487553, 
KJ487554 


this study 




Cetina 


C. dalmatina 


CETl-12 


rCETl-7 


sCETl-7 


EF605302-EF605306, 
KJ487457-KJ 487461, 
KJ487463, KJ487470, 
KJ487497-KJ 487503, 
KJ487527-KJ 487533 


this study 




Mislina (Neretva) 


C. narentana 


NERl, 3, 8-11 


rNERl-3, 5 


sNERl-3, 5, 6 


EF605315, EF605316, 
KJ487438- KJ487449, 
KJ487454-KJ 487456, 
KJ487472, KJ487485, 
KJ487507-KJ487513, 
KJ487520, KJ487521, 
KJ487539-KJ 487550 


this study 




Metkovic (Neretva) 


C. narentana 


NER5, 6, 7 


rNERl 


sNERl, 6, 10, 11 








Norin 


C. narentana 


NER2, 3, 6, 11, 12, 
14, 16, 18, 19 


rNERl, 4-6 


sNER3, 12, 13, 14 








Modro oko 


C. narentana 


NERl, 3, 11-15, 17 


rNERl, 2, 6-8 


sNERl, 5, 7-9 








Trebisnjica 


C. narentana 


NERl, 7, 12, 13 


rNERl 


sNERl-4 








Hutovo b. 


C. narentana 


NERl, 4 












Bacinska Lakes 


C. aiyrica 


BACl 


rMATl, 2 


sMATI, sKRE4 


EF605312, EF605313, 
KJ487436, KJ487437, 
KJ487451-KJ 487453, 
KJ487465, KJ587466, 
KJ487469, KJ487479- 
KJ487484, KJ487486, 
KJ487487, KJ487489- 
KJ487496, KJ487560, 
KJ487514-KJ487516, 
KJ487522, KJ487524- 
KJ487526, KJ487538, 
KJ487551, KJ487552 


this study 




Matica 


C. aiyrica 


MATl-6 


rMATl-4 


sMATI, sKREI, 4 








Krenica 


C. aiyrica 


KREl-6, MATS 


rKREl-6, rMOBI 


sKREI-4 








Prolosko b. 


C. aiyrica 


PROl-4 


rPROl-3, rKREl-3, 
rMAT4 


sPROl, sKREl 








Mostarsko blato 


C herzeQoviniensis Buj & 
Sanda, sp. nov. 


M0B1-6 


rMOBI -2, rKREl 


sMOBl, sKREl 


KJ487473-KJ 487478, 
KJ487486-KJ 487488, 
KJ487522, KJ487523 


this study 






C arachthosensis 


1 






AY191581 


[20] 








2 






AF263088 


[18] 






C. bilineata 


1-10 


rBIL2-5 


sBILl-3 


EF672359-EF672368, 
EF672415- EF672418 
EF672447-EF672449 


[22] 








11-12 






AF263090, AF263091 


[18] 










rBILI 




EF056382 


[37] 






C. calderoni 


1, 2 






AY860121, AY860122 


[64] 






C. elongata 


1-4 






EF605308, EF605310, 
EF605318, EF605320 


[21] 










rELOl 




EF056332 


[37] 










rEL02, 3 


sELOl, 2 


EF672420, EF672421 


[22] 






C. elongatoides 








AY191567 


[20] 






C. hellenica 


1, 2 






AY191582, AY191583 


[20] 
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Table 1. Cont. 





loc3lity 


Sp6Ci6S 


cyt b hsplotypss 


RAG1 

hsploty pss 


S7 hsplotypss 


3cc6ssiori nuinb6r 


rcf 6 r6 n c 6 




C. melanoleuca 




rMELl 




CF056343 


[37] 




C meridionalis 


1, 2 






AF263083, AF263084 


[18] 




C. ohridana 


1-5 






EF597224, EF597226, 
EF597227, EF597234, 
EF597240 


[65] 






6,7 


rOHRl-3 


sOHRl-3 


EF672394, AY191563, 
EF672431- EF672433, 
EF672463-EF672465 


[22] 




C. paludica 


1, 2 






AY860179, AY860180 


[64] 




C. punctilineata 


1, 2 






AY191579, AY191580 


[20] 




C. stephanidisi 


1, 2 






AY191571, AY191572 


[20] 




C. strumicae 


1,2 






DQ217372, DQ217373 


[20] 




C. taenia 


1 






AY191565 


[20] 






2 






AF263078 


[18] 








rTENl 




EF056334 


[37] 




C. tanaitica 


1, 2 






DQ217397, DQ217398 


[20] 




C. trichonica 


1, 2 






AF263085, AF263086 


[18] 




C vardarensis 


1, 2 






AY191569, AY191570 


[20] 




C. vettonica 


1, 2 






AY860182, AY860183 


[64] 



Czanandreai 1-4 rZANl, 2 sZANl, 2 EF67241 1-, EF672413, [22] 

AY191562, EF672444, 
EF672445, EF672477, 
EF672476 



doi:l 0.1 371 /journal.pone.0099833.t001 



Cobitis jadovaensis: PMF CBJAl-6, 6 spec, (morphological 
analyses: PMF CBJAl-6: 20", 49; genetic analyses: PMF CBJA 
1-6), SL 60.8-87.0 mm, Jadova R., Croatia. 

Cobitis bilimata: PMF CBZR6-32, 27 spec, (morphological 
analyses: PMF CBZR7-32: 7CT, 199; genetic analyses: PMF 
CBZR6-12, 14, 15, 16, 19, 20, 28), SL 54.0-88.0 mm, Zrmanja 
R., Croatia. 

Cobitis dalmatina: PMF CBCEl-15, PMF CBCE 18-20, PMF 
CBCE 22-26, 23 spec, (morphological analyses: PMF CBCE6-15, 

18, 19, 20, 22-26: 1 IC, 79; genetic analyses: PMF CBCEl-5, 7-9, 
12-15, 23), SL 55.8-105.3 mm, Cetina R. (in Blato village), 
Croatia. 

Cobitis mrmtam: PMF CBMIl-9, 9 spec: 40", 59 (all for 
morphological and genetic analyses), SL 56.8—88.3 mm, Neretva 
R. in MKniste village (channel Mislina), Croatia; PMF CBNO 1- 

19, 19 spec, (morphological analyses: PMF CBNOl, 3-6, 9-19: 
20-, 159; genetic analyses: PMF CBNO 1-7, 9-13, 15), SL 70.8- 
92.3 mm, Norin R., Croatia; PMF CBMOl-19, 19 spec, 
(morphological analyses: PMF CBMO 1-19: 40", 159; genetic 
analyses: PMF CBMO 1-4, 6-12, 14, 15, 18), SL 53.9-77.9 mm, 
Modro oko Lake, Croatia; PMF CBNE 5, 6, 7, 12, 13, 14, 6 spec, 
(morphological analyses: PMF CBNE6, 7, 12-14: lO", 49; genetic 
analyses: PMF CBNE5-7), SL 63.4-85.0 mm, Neretva R. in 
Metkovic, Croatia; NMP P6V 81060-81072, 81224-81239, 
81699-81704, 37 spec, (morphological analyses: NMP P6V 
81060-81072, 81224-81239: 140", 159; genetic analyses: NMP 
P6V 81699-81704), SL 50.9-74.6 mm, Trebisnjica R. (in Ravno, 
Popovo karstic field), Bosnia and Herzegovina; NMP P6V 81209- 
81223, 81928-81930, 18 .spec, (morphological analyses: NMP 
P6V 81209-81223: 4o-, 1 19; genetic analyses: NMP P6V 



81928-81930), SL 49.1-90.8 mm, Hutovo blato, Bosnia and 
Herzegovina. 

Cobitis illyma: PMF CBPR 1, 2, 4-25, 28, 29, 31, 27 spec, 
(morphological analyses: PMF CBPR7-25: lO", 189; genetic 
analyses: PMF CBPRl, 2, 4-6, 28, 29, 31), SL 52.7-98.7 mm, 
Prolosko blato Lake (Imotsko polje karstic field), Croatia; PMF 
CBBAl, 2*, 2 spec, (used only for genetic analyses), SL 58.2— 
60.4 mm, Bacinska Lakes, Croatia; PMF CBMA 1-25*, 25 spec, 
(morphological analyses: PMF CBMA 6-25: 70", 139; genetic 
analyses: PMF CBMAl-5, 10, 14-17), SL 54.0-81.9 mm, Matica 
R. (Polje Jezero karstic field), Croatia; NMP P6V 81191-81208, 
81986-81995*, 28 spec, (morphological analyses: NMP P6V 
81 191-81208: 90-, 99; genetic analyses: NMP P6V 81986-81995), 
SL 59.2-77.6 mm, Krenica Lake (Bekijsko karstic field), Bosnia 
and Herzegovina; 

Cobitis herzegovinimsis: Buj & Sanda, sp. nov: PMR VP 2950—52, 
NMP P6V 80903-80905, 81171-81190, 81959-81961, 81964, 
81967-81970, 81972-81975**, 38 spec, (morphological analyses: 
PMR VP 2950-52, NMP P6V 80903-80905, 81 171-81 190: 180", 
89; genetic analyses: NMP P6V 81959-81961, 81964, 81967- 
81970, 81972-81975), SL 51.2-73.0 mm, Mostarsko blato karstic 
field, Bosnia and Herzegovina. 

*Based on literature data, those samples were regarded as Cobitis 
narentana. 

**Based on literature data, those samples were considered as 
Cobitis narentana. 

Nomenclatural Acts 

The electronic edition of this article conforms to the requirements 
of the amended International Code of Zoological Nomenclature, 
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and hence the new names contained herein are available under 
that Code from the electronic edition of this article. This 
published work and the nomenclatural acts it contains have 
been registered in ZooBank, the online registration system for 
the ICZN. The ZooBank LSIDs (Life Science Identifiers) can be 
resolved and the associated information viewed through any 
standard web browser by appending the LSID to the prefix "http:/ / 
zoobank.org/". The LSID for this publication is: um:lsid:zoobank. 
orgipub: 261362F1-0C5A-4491-9D4A-F600075AC6BC. The 
electronic edition of this work was published in a journal with an 
ISSN, and has been archived and is available from the following 
digital repositories: PubMed Central, LOCKSS. 

Results 

The aim of the present investigation was to confirm the 
taxonomic status of all Dalmatian and Herzegovinian Cobitis 
populations. Therefore, the results were analyzed without taxo- 
nomic priors and are presented for each population. The current 
taxonomic status of each population according to the literature is 
presented in the Introduction and in Figure 1, while taxonomic 
conclusions based on our results are reported in the Discussion and 
in the other figures. 

Morphometric Characters 

Morphometric analysis was conducted separately for females 
and males since statistically significant differences were found 
between the sexes, in both measured and transformed (size- 
independent) morphometric characters (t-test; p<0.05). Mean, 
minimal and maximal values of morphometric ratios for each 
population are reported in the Table SI. 

AU size-independent morphometric characters were significant- 
ly different among both females and males from different 
populations (ANOVA; in all cases a = 0.00). However, -a post hoc 
comparison (Fisher's LSD test) revealed that some groups of 
populations are more uniform regarding their body shape than 
others. Females from Krenica, Matica (considered as C. narentana) 
and Prolosko blato (C. illyricd) possessed pectoral fins of uniform 
length, while this character was significantly different between 
those populations and the Mostarsko blato population (also 
supposedly belonging to C. narentana). Furthermore, males from 
Krenica and Matica were uniform with regard to ID, L\, IV and 
lac, whilst those measures were significandy different between 
those males and those from Mostarsko blato. 

The results of PC A analysis are shown on Figure 2. For females 
the first factor comprised 57.7% of the total variability and was 
mainly influenced by c, pan, pA, pV and pD, while the second 
factor accounted for 9.9% of the total variability and is mainly 
correlated with Ipc, IC, h and io. Similarly, in the case of males, the 
first factor encompassed 5.5.9% of the total variability and mainly 
corresponds with c, pan, pA, pV, pD and prO, while the second 
factor, which comprises 13.2°/() of the total variance, is mainly 
influenced by IC, h, io and Oh. Projection of size-independent 
measures of females revealed that the Cetina population is, based 
on standardized morphometric characters, different from all the 
other populations. Separation of populations from Prolosko blato, 
Matica and Krenica, as well as the one from Mostarsko blato, was 
also obvious. Each of these populations was, based on body shape, 
separated into its own group. Fishes from Matica R., although 
forming a special group, were included among other populations 
from the Neretva R. basin (with the exception of the Mostarsko 
blato population), which occupy the largest part of the plot. The 
plot of factor scores for males revealed a similar situation: 
specimens from Cetina, Krenica and the majority of those from 



Mostarsko blato formed separate groups, while those from the 
remaining populations were more alike. 

Results of Meristic Analysis 

In contrast with the morphometric features, the number of fin 
rays did not show any difference between females and males from 
the same population. However, differences in meristic characters 
were noticed among different populations (Table 2). The number 
of branched fm rays was constant only in pelvic fins, where 5—7 
(usually 6) rays were counted in all populations. On the other 
hand, the greatest variability was recorded in the dorsal fin rays 
number. 

Features of External Morphology 

Generally, external morphology of spined loaches from different 
rivers and lakes was quite similar. Moreover, some details of 
external morphology traditionally used to differentiate between 
species turned out to include a certain amount of intraspecific and 
intrapopulational variability (Table 3). 

Phylogeny of Dalmatian and Herzegovinian Cobitis 
Populations 

Out of 122 tissue samples used for molecular analyses we have 
obtained satisfactory cyt b sequences from 120 samples, RAGl 
sequences from 99 samples and S7 first intron from 80 samples. 
Thereafter, molecular phylogenetic analysis included 1 20 new cyt 
b sequences of length 1140 base pairs (bp), 198 nuclear RAGl 
sequences (903 bp) and 160 S7 sequences (510 bp including 
indels). For cyt b, 59 different haplotypes were found among 
Dalmatian and Herzegovinian specimens, while for RAGl and 
S7, 35 unique haplotypes were obtained for each gene. Of the 
samples analyzed, 68% were homozygous for RAGl, while only 
36% were homozygous for the S7 first intron. The greatest 
number of individuals was phased with high probabilities (= 1.00). 
Three individuals were phased with lower probabilities (0.51-0.83) 
in one or more heterozygous sites in RAGl. Inside the S7 data set, 
twelve individuals were phased with lower probabilities (0.51-0.84) 
for one or more heterozygous positions and others with high 
probabilities (>0.94). In both data sets, besides best haplotype 
guesses, all other haplotype possibihties were checked for 
ambiguous sites (with probabilities <0.94). Since ambiguous sites 
did not represent parsimony informative, nor diagnostic charac- 
ters, they turned out not to be important in phylogenetic 
reconstructions so phylogenetic analyses were based on all alleles. 
The test for base homogeneity indicated that base frequency 
distributions were always homogenous among taxa for all three 
genetic markers. Table 4 reports the phylogenetic performance of 
each gene. Neutrality tests suggested almost no deviation from 
mutation-drift equilibrium for all genes analyzed; Fu & Li's D* 
and F* statistics, as well as Tajima's D were not statistically 
significant (p>0.05 for all investigated genes and species, with 
exceptions of Fu & Li's F* and Tajima's D for cyt h in C. dalmatina, 
Tajima's D for S7 in C. narentana and Fu & Li's D* and F* for S7 in 
the population from Mostarsko blato, where 0.02>p>0.05). The 
estimated minimum number of recombination events was much 
smaller than the number of mutations in nuclear genes (5 vs. 29 in 
RAGl, 4 vs. 52 in S7 first intron). Therefore, we believe that 
phylogenetic pattern was not violated by selective forces in any of 
the investigated genetic markers. 

AU three methods of phylogenetic inference employed for cyt b 
sequences yielded trees with a similar overall topology (Figure 3). 
Samples from the investigated region, together with C. bilineata 
sequences retrieved from GenBank clustered as a monophyletic 
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Figure 2. Plot of scores for factors 1 and 2 based on size-independent measures for females (A) and males (B) from all investigated 
populations. Legend: Z = Zrmanja (C. bilineata); C = Cetina (C. dalmatina); IVle = Neretva in IVletl<ovic (C. narentana); IVIi = IVIislina (C narentana); 
N = Norin (C. narentana); IVlo = Modro ol<o (C narentana); H = Hutovo blato (C. narentana); T = Trebisnjica (C. narentana); M = Matica (C. illyrica); 
P = Prolosl<o blato (C. illyrica); K = Krenica (C. illyrica); IVIB = IVIostarsl<o blato (C. herzegoviniensis Buj & Sanda, sp. nov.). 
doi:1 0.1 371 /joumal.pone.0099833.g002 



subgroup of the so-called Adriatic phylogenetic group (smsu [21]). 
Moreover, in accordance with previous investigations [21] [22], 
the Adriatic phylogenetic group also contained sequences of C. 
elongata Heckel & Kner, 1858, C. ohridana Karaman, 1928, and C. 
zanandreai Cavicchioh, 1965. Partition of this group and the 
position of nearly all internal branches was the same in trees 
obtained by different methods, with one exception-the position of 
the haplotype from Jadova R. In previous investigations subgroups 
were denoted as "bilineata", "elongate" and "ohridana-zanan- 
dreai" clades [21] [22]. In this investigation we have analyzed 
relationships inside "bilineata" clade or subgroup. 

The abovementioned subgroup of the Adriatic phylogenetic 
group, that comprises haplotypes from the investigated area, 
consists of three evolutionary independent lineages. One, hereafter 
named as "Imotski lineage", comprised haplotypes from Matica, 



Bacinska lakes, Prolosko blato, Krenica and Mostarsko blato, 
divided into two sublineages (a division clearly visible also on the 
SP network based on cyt b sequences from these localities; see 
Figure 4). The second lineage, named hereafter "Dalmatian 
lineage", was divided into two branches. The first branch included 
all Zrmanja haplotypes, together with the Italian C. bilineata 
haplotypes, while the second branch comprised haplotypes from 
Cetina R. and Neretva R. basin (without Matica, Bacinska lakes, 
Krenica and Mostarsko blato), divided into two sublineages. The 
third lineage, named "Jadova lineage" was represented by the only 
haplotype found in the Jadova population. Based on ML and BAY 
analyses, this lineage is a sister group to the "Dalmatian lineage", 
while in the MP phylogram it represents a sister lineage to the 
"Imotski lineage". 



Table 2. Number of branched fin rays in specimens from investigated populations. 



Population/locality (species) Number of branched fin rays 





dorsal 


anal 


caudal 


pectoral 


pelvic 


Jadova (C. jadovaensis) 




6V2 (7V2) 


(14) 15 (16) 


10 (11) 


6 


Zrmanja (C. bilineata) 


6V2-7V2 


5V2 


14 


8(9) 


6-7 


Cetina {C. dalmatina) 


6V2-7'/2 


5V2 


14 (15) 


(8) 9 (10) 


6 (7) 


Metkovic (Neretva) (C narentana) 


7V2 


5V2 


14 


8(9) 


6 


Mislina (Neretva) (C. narentana) 


7V2 (8V2) 


5V2 


14 


8-9 


6 


Norin (C. narentana) 


(6V2) 7V2 


5V2 


14 


8-9 


6 


Modro oko (C. narentana) 


(6V2) 7V2 (8V2) 


5V2 


14 


8 (9) 


6 (7) 


Hutovo blato (C. narentana) 


(6V2) 7V2 


5V2 


14 


(8) 9 


6 


Trebisnjica (C. narentana) 


(6V2) 7V2 


5V2 


14 


(8) 9 


(5) 6 


Matica (C. illyrica) 


6V2 


5V2-6V2 


14 


(8) 9 


6 


Prolosko blato (C. illyrica) 


6V2 


5V2 (6V2) 


14 


(8) 9 (10) 


6 


Krenica (C illyrica) 


(5V2) 6V2 


5V2 


(13) 14 


(8) 9 


5-6 


Mostarsko blato (C. herzegoviniensis Buj & Sanda, sp. nov.) 


6V2 


5V2 (6V2) 


14 


9 (10) 


5-6 
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Table 3. The external morphology characters in investigated populations. 



Population 


Basic 


Spots on the 






suborbital 


Suborbital 


(species) 


color 


caudal fin base 


Gambetta zones 


Dorsal blotches 


spine 


spine Scale coverage 


J ad ova 


yellow 


One dark small 


Great intrapop. 


High degree of 


Hidden under 


Consist out Small, but clearly visible, 


R. (C jadovaensis) 




and oblique spot 


diversity. Various 


intrapopulational 


the skin surface. 


of two arms, scales throughout the 






in upper half of 


blotches in surface 


variability recorded 




joint in more- body surface. 






the caudal fin 


layer in Z4 {rarely 


in all populations 




less square 






base; very rarely a 


amorphous), 


varying from clearly 




angle. The 






second, oval and 


sometimes 


separated and 




longer arm is 






lighter, spot can 


merged together 


regularly shaped 




positioned 






be seen beneath 


in a stripe. Narrow, 


blotches located on 




bellow the eye 






the first one. 


paler stripe in a 


lighter or darker 




and usually ends 






Pigmentation in 


lower pigmentation 


surface, to completely 




with two 






deeper layer is 


laye 


merged ones. 




pointed 






very pale and 








branches. The 






amorphous. 








shorter arm is 
parallel with the 
anterior edge of 
the eye and 
ends with 1-3 
branches. 



Position of the 



Zrmanja 


bright 


Two clear dark 


Intrapop. diversity 


Hidden under the 


Miniature scales inserted 


R. (C bilineata) 


yellow 


oblique spots. 


noticed; squares and 


skin surface. 


into the skin. They can 






the upper one 


rectangles in Z4 




be noticed behind the 






usually darker and 


{surface layer) are 




operculum and in front 






more oblique. 


usually completely 




of the caudal fin. 






Exceptions {lower 


separated. Deeper 










spot paler and 


pigmentation layer 










vertical, or both 


usually presented by 










spots oval) 


a narrow stripe. 










extremely rare. 












Deeper layer of 












pigmentation 












absent. 








Cetina 


very pale 


Usually no spots 


Great Intrapop. 


Hidden under the 


Small and scarce scales 


R. (C dalmatina) 


yellow 


in surface layer. 


diversity. Blotches in 


skin surface 


visible on the whole 






Even when some 


surface layer in Z4 




body surface. 






kind of spots can 


various in shape, 










be noticed on the 


separated or merged 










caudal fin base. 


together. Deeper 










they are not 


layer consists of a 










darker from the 


narrow line. 










spots in the 












Gambetta zones, 












most often are 












located in the 












middle part of the 












body, irregular in 












shape and belong 












to deeper layer. 









Neretva R. in Metkovic 
(C. narentana) 



yellow 



Two spots are 
present-the upper 
one is darker and 
oblique, the lower 
one paler and 
oval (sometimes 
slightly oblique or 
irregular in shape). 
Exceptions {lower 
or both spots small 
and hard to notice, 
or upper spot 
vertical or 
irregular) are rare. 
Pigmentation in 
deeper layer 
amorphous and 
pale, often absent. 



Great intrapop. 
diversity. Surface 
pigmentation layer in 
Z4 consists of big 
rectangles or squares, 
often merged 
together. Deeper 
layer usually looks 
like a narrow, pale 
line, but sometimes is 
thicker and darker. 



Neretva R., Mislina yellow 
(C. narentana) 

Norin (C. narentana) yellow 



Largely covered 
by skin, but with 
apex sometimes 
peeking out. 



Scales pronounced, 
especially on the ventral 
parts, somewhere even 
overlapping. 



Small scales. 
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Table 3. Cont. 



Position of the 

Population Basic Spots on the suborbital Suborbital 

(species) color caudal fin base Gambetta zones Dorsal blotches spine spine Scale coverage 



Modro oko 
(C. narentana) 



yellow 



Matica R. (C illyrica) yellow 



Usually one 
dark spot 
belonging to the 
surface 
pigmentation 
layer located in 
the upper half 
of the caudal fin 
base. The 
exceptions are 
more frequent. 
Usually the spot 
is oblique; rarely 
it has a comma 
shape, or is vertical 
or oval. Very rarely 
it is small or hard to 
notice and pale, 
while exceptionally it 
cannot be noticed at 
all. Deeper layer 
developed better 
than in other 
populations, 
presented by one or 
more paler spots, 
oval or irregular in 
shape. 



Great intrapop. 
diversity. Blotches in 
the surface layer in 
Z4 rectangular, oval 
or irregular in shape; 
sometimes partly or 
completely merged 
together. Narrow, 
pale pigmentation 
line present in 
deeper layer, rarely 
thicker and/or darker. 



Mostarsko blato {C 
herzegoviniensis Buj 
& Sanda, sp. nov.) 



yellow No spots in the Intrapop. diversity 



upper layer. 
Amorphous 
pigmentation i 
deeper layer 
developed. 



noticed. Squares in 
the surface layer in 
Z4 usually partially or 
completely merged, 
sometimes forming a 
dark stripe. Deeper 
pigmentation layer 
presented by dark 
line, sometimes 
darker than blotches 
in the surface layer. 



Largely covered 
by skin, but with 
apex sometimes 
peeking out. 



Hidden under the 
skin surface. 



Scales pronounced, 
especially on the ventral 
parts, somewhere even 
overlapping. 



Trebisnjica R. 


yellow 


Greatly 


Small scales. 


C. narentana) 




uncovered. 




Hutovo blato 


yellow 


Largely covered 


Scales pronounced, 


[C. narentana) 




by skin, but with 


especially on the ventral 






apex sometimes 


parts, somewhere even 






peeking out. 


overlapping. 



Miniature scales (visible 
under magnification). 



^rolosko blato 


yellow 


Greatly 


Small scales inserted Into 


C. illyrica) 




uncovered. 


the skin. 


<renjca (C illyrica) 


yellow 


Hidden under the 


Miniature scales (visible 






skin surface. 


under magnification). 



Miniature scales (visible 
under magnification). 



doi:l 0.1 371 /journal.pone.0099833.t003 



The phylogenetic relationships recovered by RAGl haplotypes 
(Figure 5) partially diiTered from the cyt b phylogenies. Even 
though the taxon sampling was not the same (because of the fewer 
number of RAG 1 sequences of non- Adriatic species available in 
GenBank), it is surprising that the only C. taenia L., 1758 RAGl 
sequence available from GenBank clustered together with those of 
C. ohridana and C. zanandreai. Three subgroups of the Adriatic 
phylogenetic group were also recovered by RAGl, and all 
Dalmatian and Herzegovinian haplotypes were positioned inside 



one subgroup, as in the cyt h phylogenetic tree. Furthermore, 
"Imotski" and "Dalmatian lineage" were also recovered. Howev- 
er, the ' Jadova lineage" was not recovered, yet Jadova haplotypes 
clustered inside 'Tmotski lineage"; haplotypes from the Cetina 
population did not cluster inside "Dalmatian lineage" and the 
internal structuring of the lineages was not resolved during analysis 
of the RAGl gene, as it was in the cyt b phylogenies. The 
exception was C. bilineata (haplotypes from Zrmanja R. and those 
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retrieved from GenBaiik), which was also monophyletic in RAGl 
phylogenetic trees. 

In the phylogenetic reconstruction of the S7 first intron only two 
non-Adriatic species were included and they were used for rooting 
the trees. The differences between S7 phylogenetic trees (Figure 6) 
when compared to RAG 1 phylogenies are in the position of the 
Jadova haplotypes (clustered here inside the "Dalmatian lineage") 
and the monophyly of all haplotypes from Cetina R. The topology 
of the S7 phylogenetic trees was more similar to those obtained by 
phylogenetic reconstruction of cyt b haplotypes. 

A somewhat better resolution of phylogenetic relationships 
between investigated populations based on nuclear markers was 
obtained by MJ networks. In the phylogenetic network based on 
the S7 first intron (Figure 7), five clusters of haplotypes were 
separated, corresponding to the mtDNA sublineages. The 
difference is in the incorporation of haplotypes from the Mostarsko 
blato karstic field into the same cluster with the haplotypes from 
Matica, Bacinska lakes, Prolosko blato and Krenica. In the RAG 1 
haplotype network (Figure 8), two clusters corresponding to 
"Imotski" and "Dalmatian" mtDNA lineages can be recognized, 
with haplotypes from Jadova R. forming an independent lineage 
inside "Imotski" cluster and haplotypes from Zrmanja R. forming 
a monophyletic lineage inside "Dalmatian" cluster. 

It is important to mention that there was no sharing of mtDNA 
haplotypes between six sublineages recovered from the phyloge- 
netic reconstruction of the cyt b gene. On the other hand, samples 
from Mostarsko Blato, besides unique nuclear haplotypes, 
possessed RAGl and S7 haplotypes that were also found in 
individuals from Matica R., Krenica Lake and Prolosko Blato 
Lake (haplotypes rKREl and sKREl). 

Interspecific p-distances based on cyt h gene ranged from I . I to 
5.9%, while values of intraspecific p-distances were up to 0.9% 
(Table 5). Inter- and intra-specific p-distances of RAGl were lower 
than those for cyt b, as expected (Table 5). 

Bayesian Model-Testing 

Results of model averaging, namely values of Bayes factors used 
to compare estimated marginal likelihoods of possible taxonomic 
groupings, are presented in the Table 6. Based on HME and SS 
methods for marginal likelihood estimations the model A is 
preferred over all other models, whereas based on PP method the 
model C is better than the remaining models. Nevertheless, log 
BFs among models A, C and D are lower than 3, a value that is 
considered as a strong support for one model over another [58]. 

Molecular Diagnostics 

Once the taxonomic status of Dalmatian and Herzegovinian 
populations was estimated from morphological and phylogenetic 
analysis, fixed differences and shared polymorphisms among 
species were revealed in all three molecular markers investigated 
(Table 7). Diagnostic sites in cyt b were found in all species 
(Table 8) enabling molecular identification of each species. RAGl 
gene enabled molecular identification of C.jadovaensis (A at position 
444) and C. bilineata (T at 852), while based on the S7 first intron C. 
bilineata (two diagnostic sites: 487, T; 502, A) and C. dalmatina (190, 
G) could be determined. In cyt b and RAGl genes all mutations 
were substitutions. Insertions and deletions were found only in the 
S7 first intron. For the purposes of molecular diagnostics, the most 
important insertions and deletions are a single base insertion (T) at 
position 317 in C.jadovamsis and deletion of a segment comprising 
eleven nucleotides (positions 335-345) in C. narentana. Even though 
these mutations were found in most specimens of those species, 
they cannot be considered apomorphic characters because they 
were not present in all sequences. Moreover, we found individuals 
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Figure 3. ML phylogram of cyt b haplotypes showing position of and pKiylogenetic relationships among Dalmatian and 
Herzegovinian spined loaches. Numbers at nodes represent ML BS, IVIP BS and BPP values. The species delimitation (each mtDNA 
sublineage = separate species) is based on the results of this investigation, and is not completely concordant with literature data. 
doi:l 0.1 371/journal.pone.0099833.g003 



of both species that had one allele with and the other without 
mutation. 

Discussion 

Taxonomy and Distribution of Dalmatian and 

Herzegovinian Spined Loaches 

Based on phylogenetic reconstruction of mtDNA haplotypes, 
observed morphological differences and genetic distances we 
conclude that each mtDNA sublineage (with the exception of C. 
bilineata sublineages, as discussed later) represents a separate Cobitis 
species. A large number of investigations conducted on different 
fish species have found that the mutation rate of the cyt b gene 
makes it suitabk; for taxonomic studies at the level of species. The 
present investigation confirms that the phylogenetic performance 
of the cyt b gene is more adequate for such research in comparison 



with both investigated nuclear genes, because it offers a larger 
number of parsimony informative sites and exhibits a lower 
consistency index (Table 4). This conclusion was also obtained by 
Perea et al. [59] based on study of species from the Leuciscinae 
subfamily. Phylogenetic reconstructions based on nuclear genes 
are partly concordant with mtDNA phylogenies, while recorded 
discrepancies are most probably the result of incomplete lineage 
sorting. Namely, several comparative investigations have indicated 
that monophyly, or exclusivity at the level of nuclear DNA, is not a 
necessary assumption for species that have diverged more recently, 
due to the slower divergence of nuclear genes in comparison to 
mitochondrial genes [60-62]. Comparing the phylogenetic per- 
formance of two nuclear genes, we can conclude that S7 first 
intron is more appropriate for phylogenetic investigations at 
species level (due to a faster evolution rate and a greater number of 
parsimony informative characters), while RAGl is more suitable 





Figure 4. 95% parsimony network of cyt b haplotypes of spined loaches from Matica R., Bacinska Lakes, Prolosko blato, Krenica 
and Mostarsko blato. The size of ovals corresponds to haplotype frequency. Small circles are missing (unobserved) haplotypes. 
doi:1 0.1 371/journal.pone.0099833.g004 
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for higher taxonomic levels. The alternative hypothesis for 
differences in mtDNA and nuclear phylogenies, especially in the 
presence of the same nuclear haplotype in individuals belonging to 
different mtDNA lineages, would include hybridization events. 
Even though that hypothesis requires further attention, we found 
that scenario less likely. Namely, shared haplotypes were found in 
homozygote individuals of both species. Furthermore, we did not 
find a single individual with private- haplotypes of different species. 
Thereafter, we believe that the observed pattern is more consistent 
with hypothesis that the divergence of species that share some 
nuclear haplotypes happened more recently so that each still 
contains a portion of ancestral alleles. Nevertheless, a possibility of 
hybridization cannot be rejected without further investigation and 



it is also possible that both phenomena (incomplete lineage sorting 
and hybridization) can be traced in the evolutionary history of the 
Adriatic spined loaches. 

Even though we have calculated BE using three methods of 
marginal likelihood estimation, the model averaging could not be 
considered as decisive in our case due to the low BF values among 
the most probable models. Based on the results of this analysis, 
models E and F (groupings based on nuclear genes phylogenies) 
can be positively ruled out, which also speaks in favor to the 
conclusion that cyt b gene is more adequate for research of 
taxonomic relationships on the species level. On the other hand, 
positive decision regarding models A (sbc species), B (five species) 
and C (four species) cannot be made using this approach. 
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doi:1 0.1 371 /journal.pone.0099833.g006 



Nevertheless two (HME and PS) of three employed methods 
recognized model A as the most likely one, which is concordant 
with our opinion based on all other results. When PS was used for 
marginal likelUiood estimation, model C turned out as better than 
A with low BF (0.57), whereas models A and D were the same 
(BF — 0.00). In our investigation, stepping stone method for 
marginal likelihood estimation enabled the best resolution (highest 
BF factors) of the competing models. Nevertheless, even though 
this concrete case also corroborates usefulness of the recendy 
described approach using BF in species delimitation; it also 
highlights its constraints (we could not obtain positive decision 
between three most likely scenarios, PS and SS method gave 
different results) and the importance of the systematic biologists' 
critical opinion in taxonomic investigation, especially in cryptic 
species delimitation process. 



Consequendy it can be concluded that the region of Dalmatia 
and Herzegovina is in fact inhabited by six Cohith species: 

• C. jadovaensis distributed only in Jadova R. 

• C. hilineata inhabiting Zrmanja R. in Croatia. 

• C. dalmatina restricted to Cetina R. 

• C. narentana inhabiting lower parts of Neretva R. with its 
tributaries and channels (Norin and Mislina) and Modro oko 
Lake in Croatia, as well as Trebisnjica R. and the Hutovo 
blato wedand in Bosnia and Herzegovina. 

• C. illyrica distributed in Bacinska Lakes, Matica R., Prolosko 
blato and Krenica. 

• C. herzegoviniensis Buj & Sanda, sp. nov. from Mostarsko blato, a 
species described in this paper. 
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5. elongatoides 

Figure 7. Median-joining networi< of S7 first intron liaplotypes. Black circles represent median vectors. The number of mutational steps is 
displayed by the branches, when there were two or more mutations. Haplotype clusters are marked with roman numbers, corresponding to the 
numbers used to mark the mtDNA sublineages in Figure 3. The species designation is based on the results of this investigation. 
doi:1 0.1 371/journal.pone.0099833.g007 



Our results reject the assumption that C. narentana inhabits 
Matica R. [27] and reveals that its distribution is considerably 
smaller than previou.sly reported. We also found that C. illyrica is 
not re.stricted only to the Imotsko polje karstic field [15] but also 
inhabits waters of Polje Jezero karstic field in Croatia, and Bekijsko 
polje in Bosnia and Herzegovina. 

Cobitis jadovaensis is clearly separated from the other species both 
genetically and morphologically. Besides the appearance of spots 
on the caudal fin l)ase, this species differs from the others also in 
the number of branched rays in the anal, caudal and pectoral fins. 
An individual position was also supported by phylogenetic 
reconstructions based on cyt b and S7 first intron. Molecular 
diagnostic characters were found in cyt b and RAGl genes. 
Furthermore, the intraspecific p-distance of this species is much 
less than the p-cUstance between C. jadovaensis and the other 
species. 

In our study Cobitis bilineata was genetically (in all investigated 
genes) and morphologically (especially in the appearance of caudal 
fin base spots, but also in overall coloration) clearly separated from 
the remaining species. The p-distance between this and the 
remaining Dalmatian species is quite high, though phylogenetic 
analysis revealed it is more closely related to C. narentana and C. 



dalmatina. However, extraordinary deep splitting and also high 
intra.specific p-distance was recorded for C. bilineata. Namely, the 
p-distance between two sublineages comprised by this species 
(Figure 3) is ().9-1.5"/(), a level higher than intraspecific p-distances 
of all the remaining investigated species. In fact, p-distances 
recorded inside each C. bilineata sublineage resemble intraspecific 
p-distances for Adriatic spined loaches. One of the sublineages 
comprises all haplotypes from Zrmanja R. together with two 
haplotypes from the Italian Reno R, while the second sublineage 
incorporates all the remaining Italian haplotypes, including two 
haplotypes from Reno R. Therefore, these two subKneages cannot 
be considered as two different species. Their independent 
evolutionary development was probably interrupted by secondary 
contact that disabled speciation. 

The taxonomic status of spined loaches from Cetina and 
Neretva rivers has until now been uncertain due to the small 
genetic distance between them [21], [22]. However, the present 
investigation included a sufficient number of specimens to be 
assured of their independent evolutionary course and the distinct 
position of each was corroborated also by the nuclear DNA 
phylogenies. Cyt b provides diagnostic characters for both species, 
and C. dalmatina can also be determined based on the presence of 
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guanine in position 190 of tin; S7 first intron. Tiie p-distances 
between C. dalmatina and C. narentana, even thougli less than that 
between most of the species investigated, are stiU larger than the 
intraspecific p-distances. Also, fixed differences between these two 
species are present in all investigated genes, even though in lower 
numbers than among other species. Cobitis dalmatina and C. 
narentana can be distinguished based on the appearance of spots on 
the caudal fin base. Furthermore, morphometric analysis pointed 
to the separation of Cetina spined loaches from the others. 

A similar, quite small, genetic distance was recorded between C. 
illyrwa and C. herzegoviniensis Buj & Sanda, sp. nov. Nevertheless, the 
interspecific difference in cyt h is again larger than the intraspecific 
p-distances. There are also morj)hol()gical differences that separate 
these two species, especially the appearance of spots on the caudal 
fin base. The morphometric analysis implied the separation of the 
Mostarsko blato population from all the remaining ones as well. 
Furthermore, the assumption that the Mostarsko blato population 
actually represents a separate species was confirmed by the 
phylogenetic network obtained by statistical parsimony analysis of 
all cyt b sequences consisting under the "Imotski phylogenetic 
lineage". Haplotypes from Matica, Bacinska lakes, Prolosko blato 
and Krenica are separated by up to three evolutionary steps 
(mutations), while spined loaches from Mostarsko blato are 



thirteen or more steps apart (Figure 4). Although nuclear genes 
do not enable discrimination between these two species, cyt b 
provides diagnostic sites for them. Moreover, there are ten sites in 
cyt b that represent fixed differences between them. In order to 
explain the discrepancies observed between phylogenies based on 
different genes, a further investigation of the evolutionary history 
of these species is required. Nevertheless, the S7 first intron 
phylogeny offers a preliminary hypothesis: it seems that the 
sMOBl haplotype, unique to the Mostarsko blato population, is 
the basal haplotype for the "Imotski lineage" (Figure 7). The 
sKREl haplotype diverged from sMOBl and is the only 
haplotype present in both species. Based on the S7 first intron 
MJ network (Figure 7), all C. illjrka haplotypes diverged from 
sKREl, implying a possible colonization of the C. illyrica 
distribution range from the Mostarsko blato karstic field with 
indi\'iduals (■arr)'ing sKREl. Due to the slower mutation rate of 
nuclear genes in comparison with cyt b, in each species unique cyt 
b haplotypes evolved during isolation, while both species still 
contain some same nuclear haplotypes, identified as an example of 
incomplete lineage sorting. 

Among six Dalmatian and Herzegovinian Cobitis species, only C. 
bilineata has a wide distribution; all the others are endemics with 
restricted distribution ranges. Such a large number of endemic 
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species distributed in a small area, the presence of distinct species 
in neighboring localities and difTerences in the level of genetic 
diversity among species from the same area imply a complex 
evolutionary history of Cobitis in Dalmatia and Herzegovina, 
requiring further, especially population genetic, analysis. 

Diagnostic Characters for Adriatic Cobitis Species 

Morphometric characters. Until the present paper, a 
detailed comparison of morphometric characters of spined loaches 
from Dalmatia and Herzegovina had not been conducted. To 
diflferentiate between species, several authors used the ratio 
between the depth and the length of the caudal peduncle (h/ 
Ipc), the ratio between pectoral fm length and the distance between 
the base of the pectoral and pelvic fins (lp/(aV-aP)) [9], [15], head 
length in relation to its width (c/hco) and length of anal fm base in 
relation to standard length (lA/SL) [11], as well as the ratio of 
body height and standard length (H/SL), but for species outside of 
the Adriatic basin [9]. The present investigation found that none 
of these characters have a diagnostic value due to their large 
intraspecific ranges and overlap between different species (Table 
SI). In addition to small and non-representative sample sizes, 
difficulties for previous investigations lay in analyses based upon 
samples comprising both females and males where the values 
being compared largely depend on the sex ratio in the sample. 
However, the results of this investigation found no morphometric 
character to have a diagnostic value alone and therefore to be used 
to discriminate any species. On the other hand, detailed 
comparison of all standardized morphometric features pointed to 
the same taxonomic conclusions as the molecular phylogenetic 
analyses. Besides differences in body shape between some 
populations (in particular the population from Cetina R., and 
also those from Prolosko blato, Krenica and Mostarsko blato) and 
the remaining populations, PCA also led to the conclusion that the 
group comprising C. narentana populations is, based on morpho- 
metric characters, the most diverse group. Furthermore, the 
intrapopulational diversity of the Neretvanian loaches is similar to 
their interpopulational diversity, while the interpopulational 
diversity of C. illyrka is much higher than its intrapopulational 
diversity. 

Meristic characters. Even though meristic characters have 
only rarely been used for Cobitis species determination we noticed 
phylogenetic signals in the number of branched fin rays (Table 2). 
Dorsal fin rays number can enable discrimination of C. narentana 
and C. illyrica. Furthermore, based on meristic characters, C. 
jadovaensis can be separated clearly from the remaining species. 

External morphology. Among all investigated morphologi- 
cal features, only the appearance of spots on the caudal fin base 
turned out to be of diagnostic value. Although the appearance of 
Gambetta zones have been used to determine several Cobitis 
species [9], [11], [15], based on our results, the intrapopulational 
variability and similarities in the surface pigmentation layer among 
difierent species prevent the use of this character as a diagnostic 
marker. On the other hand, the deeper layer is presented with a 
narrow line in all species. Variability of the deeper layer was also 
recorded inside several populations. The same situation was 
noticed for dorsal blotches, a feature that also differs among 
specimens from the same population, with similar patterns 
recorded in the majority of populations. 

The appearance of spots on the caudal fin base is the character 
most widely used for Cobitis species determination. Karaman [26] 
mentions this feature as a diagnostic character for C. dalmatina and 
C. narentana. Since then many authors have used this feature for 
determination of several spined loaches [9], [11], [12] [14], [15]. 
Although results obtained in the present investigation are not 
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Table 6. Bayes factors (BF) for different taxonomic groupings (models) investigated in *BEAST. Bold numbers mark the highest BF 
values and the most probable model, based on each employed method (HIVlE-harmonic mean estimator, PS-path sampling, SS- 
stepping stone sampling). 
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-8.21 


1.01 




SS 


A 


B 


C 


D 


E 


F 


A 




4.94 


1.9 


2.2 


9.6 


8.34 


B 


-4.94 




-3.05 


-2.74 


4.69 


3.4 


C 


-1.9 


3.05 




0.31 


7.73 


6.44 


D 


--2.21 


2.74 


-0.31 




7.43 


6.13 


E 


-9.63 


-4.69 


-7.73 


-7.43 




-1.29 


F 


-8.34 


-3.4 


-6.44 


-6.13 


1.29 





Model A recognizes six different species; in model B C. bilineata is presented with two species (Zrmanja and Italy); in model C C. illyrica and C. herzegoviniensis Buj & 
Sanda, sp. nov. are considered as one species; model D is similar, but C. dalmatina and C. narentana are also considered as one species; model E relies on the RAG1 
phylogeny; model F incorporates S7 first intron phylogeny. 
doirl 0.1 371 /journal.pone.0099833.t006 



completely in accordance with the literature, with some intraspe- 
cific and intrapopulational diversity recorded, this feature turned 
out to be the best phenotypic diagnostic character (Table 3) 
enabling difiFerentiation of all investigated species (Figure 9). The 
characteristic spots, that enable diagnosis of several species, are 
located in the surface pigmentation layer. Pigmentation in deeper 
layer is more variable, most often amorphous and pale. It is 
present in the majority of species, absent in C. bilineata and some 
specimens of C. narentana. 



The position of the suborbital spine, namely its exposure vs. 
coverage by skin, differs among populations, even those of the 
same species. However, a constant pattern was noticed within each 
population, indicating this character's connection with ecological 
conditions at the locality, as concluded for Sabanejeivia balcanica 
(Karaman, 1922) [63]. On the other hand, we were able to reject 
the accepted view that all Cobitis species have a suborbital spine 
positioned on the skin surface [9], so that feature cannot be used 
for genus determination. 



Table 7. The number of fixed differences (regular letters) and shared polymorphisms (bold letters) among the Adriatic Cobitis 
species in cyt b/RAG1/S7 first intron. 







JAD 


BIL 


DAL 


NAR 


ILI 


HER 


JAD 




0/0/0 


1/0/0 


0/0/1 


0/0/1 


0/0/0 


BIL 


46/6/4 




0/1/0 


1/0/0 


1/0/0 


0/0/0 


DAL 


34/3/2 


34/2/6 




0/0/1 


2/0/0 


0/0/0 


NAR 


36/6/2 


38/1/2 


7/1/3 




4/0/0 


0/0/0 


ILI 


36/2/0 


53/4/4 


43/2/2 


44/4/2 




0/3/0 


HER 


46/2/2 


56/4/6 


49/2/4 


48/4/4 


10/0/0 





Abbreviations: }AD = C. jadovaensis; BIL = C bilineata; DAL = C. dalmatina; NAR = C narentana; ILI = C illyrica; HER = C. herzegoviniensis Buj & Sanda, sp. nov. 
doirl 0.1 371/journal.pone.0099833.t007 
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Table 8. Diagnostic characters in the cyt b gene for the 
Adriatic Cobitis species. 



species 


position 


base 


C jadovaensis 


102 


T (vs. C, A) 




444 


C (vs. T) 




588 


T (vs. C) 




642 


T (vs. C) 




960 


C (vs. T) 




982 


T (vs. C) 




1047 


C (vs. T) 


C. bilineata 


201 


C (vs. T) 




213 


C (vs. T) 




265 


C (vs. T) 




276 


T (vs. C) 




360 


C (vs. T) 




477 


T (vs. C) 




478 


T (vs. C) 




501 


T (vs. A, G) 




546 


T (vs. C) 




600 


T (vs. C) 




631 


C (vs. T) 




675 


G (vs. A) 




924 


C (vs. T) 




987 


T (vs. C) 




1023 


A (vs. G) 


C. dalmatina 


447 


T (vs. C) 


C. narentana 


52 


T (vs. C) 




219 


T (vs. G, A) 




267 


C (vs. T) 




580 


C (vs. A) 




897 


A (vs. G) 


C. itlyrica 


352 


G (vs. A) 




1005 


G (vs. A) 


C. herzegoviniensis Buj & Sanda, sp. nov. 


132 


G (vs. A) 




318 


C (vs. T, A) 




327 


C (vs. T) 




660 


C (vs. T) 




954 


G, C (vs. A) 



doi:l 0.1 371 /journal.pone.0099833.t008 



Molecular diagnostic characters. For the purpose of 
molecular diagnostics and barcoding, we propose utilization of 
the gene for cytochrome b that provides diagnostic characters for 
all Adriatic Cobitis species (summarized in Table 8). 

Contribution to Fish Taxonomy 

In contributing to the taxonomy and determination of Cobitis 
species, based on the results obtained here, we bring a description 
of a new species from Mostarsko blato karstic field and a diagnosis 
for other species from the investigated area, as well as a key for 
their determination. 

Cobitis herzegoviniensis Buj & Sanda, sp. nov. urnilsid:- 
zoobank.org:act: 24564E1D-1D17^D9A-AD41-B9D030D06CB3 
(Figures 10 and 11). 



Holotype. NMP P6V 80904; CT; 55.65 mm SL, the Tistica 
River in the Mostarsko blato karstic field, Bosnia and Herzego- 
vina, N43°22'41.3" E17°19'98.3" elevation 227 m (leg. R. Sanda, 
A. Perdices, S. Perea, I. Bogut; 23.07.2004). 

Paratypes. NMP P6V 80903, 80905, PMR VP 2950-52; 
data as for holotype; NMP P6V 81173, 75-77, 79; artificial 
channel in Mostarsko blato karstic field, Bosnia and Herzegovina, 
N43°19.360' E17°43.099' elevation 222 m (leg. R. Sanda, J. 
Kohout, A. Sediva, L Bogut, 13.07.2006); SO" and 59; 44.2- 
63.96 mm ST (Figure 12F). 

Type locality. Tistica River in Mostarsko blato, Bosnia and 
Herzegovina. 

ADDITIONAT MATERMT: NMP P6V 81171, 72, 74, 78, 
81180-90; 120- and 39; 51.47-76 mm ST; the same data as for 
paratypes NMP P6V 81173, 75-77, 79. 

Etymology. The species name represents the name of the 
province where Mostarsko blato is located; an adjective. 

Diagnosis. Cobitis herzegoviniensis Buj & Sanda, sp. nov. is 
distinguished from other species of the genus Cobitis by the 
following combination of characters: a single Canestrini scale on 
the pectoral fins of males, triangular in shape with roundish outer 
border; usually b / 2 branched anal fin rays, 6/2 branched dorsal 
fin rays, 1 4 branched caudal fin rays; no characteristic spots on the 
caudal fin base, but pale deeper pigmentation layer, irregular in 
shape, can sometimes be noticed; scales on the body very small, 
visible only under magnification. 

Description. D II(III)+6V2, A II-III+5V2(6 Vj), P 9(10), V 
1+5-6, C I+14+I. Body small, larger in females (TT 60-76 mm in 
females, 51-66 mm in males; ST 52—67 mm in females, 44- 
56 mm in males), very elongated. Body depth uniform throughout 
the whole body length; maximum body depth (in front of dorsal fin 
origin) 11.4-16.3% of ST. Caudal peduncle depth 47.2-71.5% of 
its length in females; 40.9-69% in males. Head very elongated, 
comprising 19.6-23% of ST. Very small (visible only under 
magnification), separated scales distributed throughout the body 
surface. Suborbital spine located below eyes and completely 
covered by a skin fold. Three pairs of barbels located around the 
mouth, variable in size. Pelvic fins shorter than the distance 
between pelvic fin base and anal aperture. Straight margins of 
dorsal and anal fins. 

Males different from females in having longer pectoral fins (in 
males they occupy 65-78% of distance between pectoral and 
pelvic fins origins, while in females they occupy 28-45%), longer 
pelvic fins (10.8-15.6% of ST in males vs. 9.8-10.8% of ST in 
females), longer caudal fin (14.1-16.7% of ST in males vs. 12.5— 
15.4% of ST in females), longer head (19.6-23.0% of ST in males 
vs. 19.9-21.8% of ST in females) and the presence of a single 
Canestrini scale on a dorsal surface of pectoral fins. All the 
remaining body dimensions are larger in females than in males. 

Coloration (preserved specimens) . Varying from distincdy 
yellow with dark black spots and blotches to somewhat paler 
overall coloration. The ventral parts yellow without pigmentation, 
or with few spots on fins. Small spots located also on the head, 
together with a stripe that usually ends in front of the eye, and 
rarely crosses the eye. Spots are often larger on the operculum and 
on the dorsal part of the head. Spots on the caudal fin base are not 
clearly visible and usually cannot be distinguished from the 
remaining body spots; even in specimens where spots (one or two) 
are noticed, they are small and hardly visible. On the body sides 
four Gambetta zones (Z1-Z4) can be recognized. Zl is most often 
narrow and not pronounced, and is composed of very small, dense 
spots. In some specimens those spots are miniature and pale, even 
to the extent that they are hardly visible. Z2 consists of rectangular 
blotches that can be separated, but densely spaced, or merged into 
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Figure 9. Pictures of the caudal parts of specimens from investigated Cobitis species. Legend: A = C. jadovaensis; B = C. bilineata; C = C. 
dalmatina; D = C. narentana; E = C. illyrica; F = C herzegoviniensis Buj & Sanda, sp. nov. Characteristic spots on the base of the caudal fin can be seen. 
doi:1 0.1 371 /journal.pone.0099833.g009 



the single line. Z3 most often contains small pale spots. Often there 
are no spots in the posterior part of the Z3 and rarely Z3 is 
composed of only a few spots behind the head. Spots in Z3 are 
often located in a single line, rarely in several lines. In the majority 
of specimens Z4 is fiUed with square or rectangle blotches 
(presenting the surface pigmentation layer) on a darker line 
(located in the deeper pigmentation layer) that is usually more 
pronounced than blotches, while in posterior parts the blotches are 
interconnected into a single line. Sometimes, blotches are very 
pale and only a dark line is clearly visible. Dorsally, a single row of 
blotches is located, often irregular in shape and only sometimes 
squared. Dorsal blotches are sometimes separated, but located 
close to one another, and sometimes joined. Sexual dimorphism 
was not noticed in the coloration. 

Comparative remarks. Cobitis herzegoviniensis Buj & Sanda, 
sp. nov. can be distinguished from C. illyrica by the appearance of 
spots on the caudal fin base (no pronounced spots in C. 
herzegoviniensis Buj & Sanda, sp. nov. vs. one spot in the surface 
pigmentation layer of the caudal fin base in C. illyrica (exceptionally 
those spots are hardly visible). Cobitis herzegoviniensis Buj & Sanda, 
sp. nov. differs from C. narentana in the number of branched rays in 



the dorsal fin (6 in C. 



Buj & Sanda, sp. nov. vs. 



nov. vs. two clearly visible spots in C. narentana). Cobitis 
herzegoviniensis Buj & Sanda, sp. nov. is different from C. dalmatina 
in the number of branched rays in the pelvic (5-6 in C. 
herzegoviniensis Buj & Sanda, sp. nov. vs. 6 (7) in C. dalmatina) and 
dorsal fin (6' in C. herzegoviniensis Buj & Sanda, sp. nov. vs. 6V2- 



usually 7 /2 in C narentana), the appearance of spots on the caudal 
fin base (no pronounced spots in C. herzegoviniensis Buj & Sanda, sp. 



7 /2 in C. dalmatina). Cobitis herzegoviniensis Buj & Sanda, sp. nov. 
can also be distinguished from C. dalmatina by the appearance of 
scales (miniature, visible only under magnification in C. herzegovi- 
niensis Buj & Sanda, sp. nov. vs. small, but clearly visible in C. 
dalmatina). Cobitis herzegoviniensis Buj & Sanda, sp. nov. differs from 
C. bilineata in the appearance of spots on the caudal fin base (no 
pronounced spots in C. herzegoviniensis Buj & Sanda, sp. nov. vs. two 
clearly visible spots in C. bilineata). Differences between C. 
herzegoviniensis Buj & Sanda, sp. nov. and C. jadovaensis concern 
the appearance of spots on the caudal fin base (no pronounced 
spots in C. herzegoviniensis Buj & Sanda, sp. nov. vs. one spot in C. 
jadovaensis) and the number of branched rays in the dorsal (6 ' / j vs. 
7'/2), anal (usually 5V2 vs. usually 6V2), caudal (14 vs. usually 15) 
and pectoral fins (usually 9 vs. usually 10). 

Distribution. C. herzegoviniensis Buj & Sanda, sp. nov. is 
endemic to Mostarsko blato karstic field in Bosnia and 
Herzegovina. 

cobitis jadovaensis Mustafic & Mrakovcic, 2004 (Figure 
12A). 




Figure 10. Cobitis herzegoviniensis Buj & Sanda, sp. nov., a newly described species from Mostarsko blato; holotype. 

doi:1 0.1 371/journal.pone.0099833.g01 0 
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Figure 1 1 . Photograph of a live specimen of C. herzegoviniensis Buj & Sanda, sp. nov. 

doi:1 0.1 371 /joumal.pone.0099833.g01 1 




Figures 12. Photographs of representative individuals for each 
species. Legend: h = C. jadovaensis; B = C. bilineata; C = C. dalmatina; 
D = C. narentana; E = C. illyrica; F = C. herzegoviniensis, paratype PIVIR VP2950. 
doi:1 0.1 371/journal.pone.0099833.g01 2 



Diagnosis. Cohitis jadovamsis is distinguished from other 
species of the genus by the following combination of characters: 
a single triangular with roundish outer border Canestrini scale on 
the pectoral fins of males; usually S'/j branched anal fin rays, 7V2 
branched dorsal fm rays, usually 15 branched caudal fin rays, 
usually 10 branched pectoral fin rays; one small dark oblique spot 
located in the upper half of the caudal fin base, pigmentation in 
deeper layer not pronounced and irregular in shape; whole body 
covered with small, yet clearly visible, scales. 

Comparative remarks. Cobitis jadovaensis differs from C. 
bilineata, C. dalmatina, C. narentana, C. illyrica and C. herzegoviniensis 
Buj & Sanda, sp. nov. by having usually 6V2 (sometimes 7V2) 
branched anal fm rays (vs. usually 5V2); usually 15 branched 
caudal fin rays (vs. usually 14); usually 10 (vs. usually 8 or 9) 
branched pectoral fin rays. Cobitis jadovaensis can be distinguished 
from C. bilineata, C. dalmatina, C. narentana and C. herzegoviniensis Buj 
& Sanda, sp. nov. in the appearance of spots on the caudal fin base 
(usually one spot in C. jadovaensis vs. two spots or no pronounced 
spots). 

Distribution. C. jadovaensis is endemic to Jadova River in 
Croatia. 

cobitis bilineata Canestrini, 1865 (Figure 12B). 

Diagnosis. Cohitis bilineata is distinguished from the other 
species of the genus by the following combination of characters: a 
single triangular with roundish outer border Canestrini scale on 
the pectoral fins of males; 5/2 branched anal fin rays, 14 
branched caudal fin rays, usually 8 branched pectoral fm rays; two 
clear dark oblique spots on the caudal fin base, the upper one 
usually darker and more oblique; miniature scales inserted into the 
skin. 

Comparative remarks. Cobitis bilineata is difiFerent from C. 
jadovaensis, C. dabnatina, C. illyrica and C. herzegoviniensis Buj & Sanda, 
sp. nov. in the appearance of the spots on the caudal fin base (two 
spots in C. bilineata vs. usually one or no pronounced spots). Cobitis 
bilineata can be distinguished from C. narentana by the appearance 
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of the lower spot on the caudal fin base (usually obhque in C. 
bilineata vs. usually oval in C. narentam); overall coloration (bright 
yellow with clear dark spots in Gambetta zones in C. bilineata vs. 
paler body color and paler spots in C. narentana); and the 
appearance of spots in Z4 (usually clearly separated, big squares 
and rectangles in the surface pigmentation layer of the Z4 in C. 
bilineata vs. often partly or completely merged squares and 
rectangles, sometimes oval in shape in C. narentana). 

Distribution in Croatia. Zrmanja R. The distribution range 
of C. bilineata outside of Croatia was not investigated. 

cobitis dalmatina Karaman, 1928 (Figure 12C). 

Diagnosis. Cobitis dalmatina is distinguished from the other 
species of the genus by the following combination of characters: a 
single triangular with roundish outer border Canestrini scale on 
the pectoral fins of males; usually 5V2 branched anal fin rays, 
usually 14 branched caudal fin rays, usually 9 branched pectoral 
fin rays; no pronounced spots on the caudal fin base (no spots in 
the surface layer, but deeper pigmentation layer usually present); 
small scales on the whole body surface, but scarce. 

Comparative remarks. Cobitis dalmatina can be distinguished 
from C. jadovaensis, C. bilineata, C. narentana and C. illyrica by the 
appearance of spots on the caudal fin base (no pronounced spots in 
C. dalmatina vs. 1 or 2 spots). Cobitis dalmatina dilfers from C. 
herzegoviniensis Buj & Sanda, sp. nov. in having 6-7 branched pelvic 
fin rays (vs. 5-6) and in having small, but clearly visible, scales on 
the body surface (vs. miniature scales visible only under 
magnification). 

Distribution. Cobitis dalmatina is endemic to Cetina R. in 
Croatia. 

cobitis narentana Karaman, 1828 (Figure 12D). 

Diagnosis. Cobitis narentana is different from the remaining 
species of the genus by the following combination of characters: a 
single Canestrini scale on the pectoral fins of males, triangular in 
shape with roundish outer border; 5/2 branched anal fin rays, 
usually 7V2 branched dorsal fin rays, 14 branched caudal fin rays, 
8-9 branched pectoral fin rays; two spots on the caudal fin base- 
the upper one darker and oblique, the lower spot paler and usually 
oval; scales small or pronounced, but always visible. 

Comparative remarks. Cobitis narentana differs from C. 
jadovaensis, C. dalmatina, C. illyrica and C. herzegoviniensis Buj & 
Sanda, sp. nov. in having usually 2 spots on the caudal fin base (vs. 
one spot or no clearly visible spots). It can be distinguished from C. 
Ulrica and C. herzegoviniensis Buj & Sanda, sp. nov. by the number 
of branched dorsal fin rays (usually 7V2 in C. narentana vs. usually 
6V2). Cobitis narentana is different from C. bilineata in the 
appearance of the lower spot on the caudal fin base (usually oval 
in C. narentana vs. oblique in C. hilineata); overall coloration (paler 
body color and paler spots in C. narentana vs. bright yellow body 
color with clear dark spots in Gambetta zones in C. bilineata); and 
the appearance of spots in Z4 (often partly or completely merged 
squares and rectangles, sometimes oval in shape in the surface 
pigmentation layer of the Z4 in C. narentana vs. usually clearly 
separated, big squares and rectangles in C. bilineata). 

Distribution. C. narentana is distributed in the lower parts of 
Neretva R. with its tributaries and channels and Modro oko Lake 
in Croatia, as well as in Trebisnjica R. and Hutovo blato wetiand 
in Bosnia and Hczego\ ina. 

cobitis illyrica Frcyhof & Stelbrink, 2007 (Figure 12E). 

Diagnosis. Cobitis iltfrica can be distinguished from the other 
species of the genus by the following combination of characters: a 
single triangular with roundish outer border Canestrini scale on 
the pectoral fins of males; usually 5V2 branched anal fin rays, 
usually 6V2 branched dorsal fin rays, usually 14 branched caudal 
fin rays, usually 9 branched pectoral fin rays; usually one oblique 



spot in the upper half of the caudal fin base; small or miniature 
scales inserted into the skin. 

Comparative remarks. Cobitis illyrica differs from C. 

bilineata, C. dalmatina, C. narentana and C. herzegoviniensis Buj & 
Sanda, sp. nov. in having usually one spot on the caudal fin base 
(vs. two spots or no spots). It is different from C. narentana in the 
number of branched dorsal fin rays (usually 6V2 vs. usually 7'/2). 
Cobitis illyrica is different from C. jadovaensis in the number of 
branched dorsal fin rays (5V2-6V2 in C. illyrica vs. 7V2 in C. 
jadovaensis); the number of branched caudal fin rays (usually 14 in 
C. illyrica vs. usually 15 in C. jadovaensis); the number of branched 
pectoral fin rays (usually 9 in C. illyrica vs. usually 10 in C. 
jadovaensis). 

Distribution. Cobitis illyrica inhabits Matica R., Prolosko 
blato and Bacinska Lakes in Croatia, as well as Krenica Lake in 
Bosnia and Herzegovina. 

Determination Key to Cobitis Species in Dalmatia and 
Herzegovina 

la The number of branched fin rays in the anal fin 6V2 
(exceptionally 7/2) and in the caudal fin 1 5 (exceptionally 1 4 and 
16). The whole body covered with small, yet clearly visible scales. In 
the surface pigmentation layer of the caudal fin base usually one 
small, oblique, dark spot located in the upper half of the body sides. 
Cobitis jadovaensis 

lb In the anal fin usually 5V2 (exceptionally 6V2) and in the 



caudal fin 14 branched rays 2 

2 a Two dark spots in the surface pigmentation 

layer of the caudal fin base 3 

2b One dark spot or no dark spots in the surface 

pigmentation layer of the caudal fin base 4 



3a On the caudal fin base clearly visible two dark spots (located 
in the surface pigmentation layer), most often oblique in shape, the 
upper one being larger. In Z4 usually clearly separated big squares 
and rectangles. In dorsal fin 6V2-7V2 branched rays. 
C. bilineata 

3b Usually two dark spots in the surface pigmentation layer of 
the caudal fin base-the upper one darker and oblique, the lower 
one being lighter and roundish, even though exceptions are 
possible. Squares and rectangles in Z4 often partly or completely 
merged. In dorsal fin usually 7V2, exceptionally 6V2 or 8V2 
branched rays C. narentana 

4a On the caudal fin base only an upper dark spot is usually 
visible (belonging to the surface pigmentation layer); exceptionally 
anoth(-r, 1()W(t and lighter, spot (belonging to the deeper 
pigmentation layer present. In dorsal fin 6V2, exceptionally 5V2 
branched rays C. illyrica 

4b No pronounced spots on the caudal fin base, rarely 
small hardly visible spots irregular in shape and belonging to the 
deeper pigmentation layer can be noticed. The number of 
branched dorsal fin rays 6V2-7V2 5 

5a No dark spots on the caudal fin base. The number of 
branched dorsal fin rays 6V2-7V2. Small scales on the whole 
body surface C. dalmatina 

5b No spots on the caudal fin base, or, rarely, small and hardly 
visible spots present, located in the deeper layer. The number of 
branched dorsal fin rays 6'/2. ScEiles very small, visible only under 
magnification C. herzegoviniensis Buj & Sanda, sp. nov. 

Supporting Information 

Table SI Mean, minimal and maximal values of 
morphometric ratios of investigated Cobitis species. 
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